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The conservation of the number of particles within the QRPA plays an important role in the 
evaluation muon capture rates in all light nuclei with A < 30. The violation of the CVC by the 
Coulomb field in this mass region is of minor importance, but this effect could be quite relevant 
for medium and heavy nuclei studied previously. The extreme sensitivity of the muon capture 
rates on the pp coupling strength in nuclei with large neutron excess when described within the 
QRPA is pointed out. We reckon that the comparison between theory and data for the inclusive 
muon capture is not a fully satisfactory test on the nuclear model that is used. The exclusive muon 
transitions are much more robust for such a purpose. 
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1. Introduction 

Among different semileptonic processes, the muon capture is one of the weak observables 
that, together with the /3 -decay, has available a fruitful set of experimental data that were collected 
in the last fifty years. Several works were focused to establish the universal V — A character of 
nuclear muon capture, the role of induced currents, second-class currents, and nonexistence of 
V +A interactions. It is known that the experimental value of the induced pseudoscalar coupling 
g P is the least known of the four constants (g v ,gA,gM,gp) defining the weak nucleon current. 

Its size is dictated by chiral symmetry arguments, and its measurement represents an important 
test of quantum chromodynamics at low energies [JTJ . During the past two decades a large body 
of new data relevant to the coupling g P has been accumulated from measurements of radiative and 
non radiative muon capture on targets ranging from 3 He to complex nuclei. Only transitions to 
unnatural parity states depend on g P , as can be seen from Eq. ([Q|). A summary of references on 
these issues are cited in review papers Ref. [Q [5| 0. 

Simultaneously, the muon capture processes have been used to scrutinize the nuclear structure 
models, since they provide a testing ground for wave functions and, indirectly, for the interactions 
that generate them. Being the momentum transfer of the order of the muon mass m^ = 105.6 MeV, 
the phase space and the nuclear response favor lower nuclear excitation energies, and thus the tran- 
sitions to nuclear states in the giant resonance region are the dominate ones. We will cite only a few 
of them. Most of these works were done within the shell model (SM) framework [Q, ||, |6[ ^ pj. Sev- 
eral studies were performed by employing the random phase approximation (RPA) [|7| ^ |10|, O], 
In the last work, where the total muon capture rates for a large number of nuclei with 6 < Z < 94 
have been evaluated, the authors claimed that an important benchmark was obtained by introducing 
the pairing correlations. They have done this ad-hoc by multiplying the one-body transition ma- 
trix elements by the BCS occupation probabilities. However, we know that the quasiparticle RPA 
(QRPA) formalism is a full self-consistent procedure to describe consistently both i) short-range 
particle-particle (pp) pairing correlations, and ii) long-range particle-hole (ph), correlations han- 
dled with RPA. Quite recently, the relativistic QRPA (RQRPA) was applied in the calculation of 
total muon capture rates on a large set of nuclei from 12 C to 244 Pu, for which experimental values 
are available [p"3|]. 

In the present work we do a systematic study of the muon capture rates of nuclei with 12 < 
A < 56 masses ( 12 C, 20 Ne, 24 Mg, 28 Si, 40 Ar, 52 Cr, 54 Cr, 56 Fe) within the number projected QRPA 
(PQRPA). The motivation for this investigation comes from the successful description of weak 
observables in the triad { 12 B, 12 C, 12 N} within this model [|TJ, 14]. There, it was shown that the 



projection procedure played an essential role in properly accounting for the configuration mixing in 
the ground state wave function of 12 N. The employment of PQRPA for the inclusive 12 C ( v e , e~ ) 12 N 
cross section, instead of the continuum RPA (CRPA) used by the LSND collaboration in the anal- 
ysis of Vj, — > v e oscillations of the 1993-1995 data sample, leads to an increased oscillation prob- 
ability tJ15H. The charge-exchange PQRPA, derived from the time-dependent variational principle, 



was used to study the two-neutrino /3/3 -decay amplitude ^2v in 76 Ge [16]. In that work, the pro- 
jection procedure was less important and the QRPA and PQRPA yield qualitatively similar results 
for ^#2v- The PQRPA was recently used to calculate the 56 Fe(v e ,e~) 56 Co cross section [17]. 



We will also give a glance on the violation of the CVC by the Coulomb field, which was 
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worked out recently [14], and appears in the first operator ( 2 .4 ) for natural parity states i . This effect 
is expected to be tiny for the nuclei studied here, since ABcoul is relatively small in comparison with 
m M ; it goes from 3.8 MeV in 12 C to 9.8 MeV in 56 Fe. 

2. /i -capture rates formalism 

When negative muons pass through matter, they can be captured into high-lying atomic or- 
bitals. From there they then quickly cascade down into the \S orbit with binding energy Eg, 
where two competing processes occur: one is ordinary decay n~ — > e~ + v^ + v e with character- 
istic free lifetime 2.197 x 10 6 sec, and the other is (weak) capture by the nucleus }i~ + (Z,A) — > 
(Z — 1 , A) + v^. The latter, naively expected to scale with Z, is drastically enhanced by an additional 
factor of Z 3 , originating from the square of the atomic wave function <p\s evaluated at the origin 
[2]. Thus, its rate is roughly proportional to Z 4 and dominates decay at large Z. This dominance is 
however significantly diminished by the gradual decrease of the effective-charge correction factor 



®{Z) p, 18] 



Then the muon capture rate from the ground state in the initial nucleus (Z,A) to the state J^ in 
the final nucleus (Z—l,A) reads 

A(J«) = §|0is| 2 ^(Z)^ A (£v,J*), (2-1) 

where 

E v = k = m M - (M„ - M p ) -Eg- Oj* (2.2) 

is the neutrino energy, and 

^ A (£ v ,J^)=47rG 2 [|(J^||O 0J (£ v )-Ooj(£v)||O + )| 2 + 2|(J^||O_ 1 j(£ v )||O + )| 2 ], (2.3) 

is the transition probability, being the Fermi coupling constant G = (3.04545 ± 0.00006) x 10~ 12 
natural units. The nuclear operators are: 



<~>0J - O J = gy ^»j , 



E v 

y.\J 



O-ij = -{gK + Fw)^_ij +gv^T u , (2.4) 

for natural parity states (n = (—) J , i.e., J K = + , 1 ,2 + , 3 , • • •), and 

O j - O J = gA-^f + C?A + g~A~ &>) -^OJ I 

O-u = -(gA + g^^j-gy^Zj, (2-5) 

for unnatural parity states (n = (— ) /+1 , i.e., J n = 0~, l + ,2~,3 + , • • •). The elementary operators 



are: 



jq = js{p)Ys{i) , ^j=M- 1 £/ J -'- 1 F LJm A(p)[7 L (f)®V] J , 

L>0 

Jtj =M- 1 jj(p)F J (r)(a.V) , ^ = £ i J - L - l F Um j L (p) [F L (r) ®a]j (2.6) 

L>0 



'When the consequences of the CVC are not considered, as in Ref. [J3|], the factor (m^ — A£c ou i —E^)/E v in this 
relation goes to unity. 



Muon capture within the PQRPA Danilo Sande Santos 



where F|_j m = (— ) (1, — mSm\ LO) is a Clebsch-Gordon coefficient, p = |k|r, and the superscripts 



R, and / in ( [2 .4] ) and ( |2.5| ) stand for real and imaginary pieces of the operators ( |2.6| ). Moreover, 

Ai'coul = 4^- - 1-45ZA" 1/3 MeV, E% = (eZ) 2 ^- =■ 2.66 x 10~ 5 Z 2 mu, (2.7) 

5/v 2 

and 



\"^T7' Kw i.S'v + ^Mj^rr; g P — ^p^TTj (2-8) 



v — _ /■ I \ v — _ v 
V 2M' ? *-« v + W2M : 5p_5p 2M 

with g v , g A , g M , and g P , being the effective vector, axial-vector, weak magnetism, and pseudo-scalar 

coupling constants, respectively. We adopt 

2Mm u 
g v = l, g A = 1.135, g M = 3.70, g P = g A j- 2 \ =6.7, (2.9) 

yt 2 + m2 

where the value for g P comes from the PCAC, pion-pole dominance and the Goldberger-Trieman 



relation [19], and for g A we use the same value as in Ref. [13]. 



From Eqs. (A6) and (A7) in Ref. [ |14] ] one sees that g P is contained in axial-vector pieces of 
both operators O^j (temporal), and Ooj (spacial). They contribute destructively, being dominant 
the second one. In Ref. [|13|] g P appears only in the temporal operator. However, after making use 
of the energy conservation condition (jl2|), i.e., K = m^ +k@ (k® = — Wjj) one ends up with the 
same result for Obj — Ooj- 

The + «->■ _ transitions are determined by two nuclear matrix elements only: ^#q and ^q , 



as can be seen from the first relation in (2.5). As such they are the most appropriate to extract 
the magnitude of g P from the muon capture experiments. In fact, studies of the 16 O(0j") — > 16 
N(0j~) transition within large-basis SM calculations have yielded values of g P = 6 — 9 [g], and 
g P = 7.5 ± 0.5 [§] that are consistent with the estimate (|^) as well as with theoretical prediction 
g P = 8.2 from chiral symmetry arguments [p. More recently, Gorringe [Q] reported from the SM 
study of muon capture 40 Ca(0 ] F ) — > 40 K(O^) have extracted from the experimental result of A the 
values g P = 14.3+] ;|, and g P = 10.3+^. 

3. Numerical results 

For the set of nuclei discussed here we have adopted the single-particle energies (s.p.e.) from 
the self-consistent calculation performed by Marketin et al. Jl3| ] within the relativistic Hartree- 
Bogoliubov model (RHB), using effective Lagrangians with density-dependent meson-nucleon 
couplings and DD-ME2 parametrization. The residual interaction is approximated by the 5-force 
(in MeV fm 3 ) 

V = -4n(v s P s + v t P,)S(r), 

with singlet (v s ), and triplet (v t ) coupling constants different ph, pp, and pairing channels. The 
proton and neutron pairing parameters v^ air {p) and v p s a,r {n), used in solving the BCS and PBCS 
equations, were determined from the experimental data by the adjusting procedure described in 
Ref. [20]. For the parameters in the ph channel we employe the values vf = 27 and vf = 64, 
which were obtained in a systematic study of the GT resonances [21 ]. The ratio t = vf p /v p s mr was 



considered as free parameter within the pp channel. It was found [22] that the muon capture, just 
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Figure 1: (Color online) Ratios of theoretical to experimental inclusive muon capture rates for different 
nuclear models, as function of the mass number A. The present QRPA and PQRPA results, as well as 
the RQRPA calculation Q were done with g A = 1.135, while in the RPA+BCS model [[Ti}| was used the 
unquenched value g A = 1.26 for all multipole operators, except for the GT ones where it was reduced to 

J?A~1. 



like the j3/3-decay, probes the final leg of a /3/3 -transition and as such strongly depends on the 
strength of the pp interaction. Even worse, the QRPA model collapses as whole in the physical 
region of t pl^, 21, 23]. Yet, the distinction between the initial and final legs in the j3j3 -decay only 
makes sense in nuclei that possess an appreciable neutron excess, which doesn't happen in nuclei 
under discussion where N = Z. Moreover, the results of the PQRPA calculations in 12 C, displayed 
in Fig. 5 of Ref. [|4|] suggest that the choice t = could be appropriate for the description of N = Z 
nuclei. Therefore, this value of the pp coupling strength is adopted here. 

Ratios of theoretical to experimental inclusive muon capture rates for different nuclear models 
are exhibited in Fig. ||. It is self evident that the number projection plays an important role in 
light nuclei with A ^ 30, making that the PQRPA agrees better with data than the plain QRPA. On 
the other hand it is difficult to judge whether our estimates are better or worse than the previous 



ones nil, 13]. 



We have found that the consequences of the violation of the CVC by the Coulomb poten- 
tial [Q] for the nuclei considered here is very tiny. In fact, the major effect appears in 56 Fe, where 
the total muon capture is reduced from A = 4260 x 10 3 s~ ! to A = 4056 x 10 3 s _1 . 

In the case of 12 C we have at our disposal also the experimental data for exclusive muon 
capture rates to bound excited states J^ = 1^,2^,2^, and lj~ in 12 B |§, 26]. They have been 



discussed previously in the framework of the PQRPA [ |12] [24] ], but for the sake of completeness 
we show them again in Table [j] The most relevant to highlight in this table is that, while both 
PQRPA calculations of the inclusive muon capture rates agree fairly well with the experiment, 
the corresponding exclusive reactions are very different in the two calculations. In other words, the 
agreement between theory and data for the inclusive muon capture does not guarantee the goodness 
of the model that is used. 
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Table 1: Energies (in units of MeV) and exclusive muon capture rates (in units of 10 3 s ) for the bound 
excited states in 12 B. Besides the present PQRPA result, we also show a previous one [|l2||, as well as those 
evaluated within the SM [§], and the RPA ||[[l0||. 



Model 


J* 


n 


2+ 


2r 


ir 


J^inc 


PQRPA 


E 


0.00 


0.43 


6.33 


6.83 






A 


8.80 


0.20 


0.60 


0.85 


37 


PQRPA [12] 


E 


0.00 


0.50 


2.82 


3.31 






A 


6.50 


0.16 


0.18 


0.51 


40 


SM[§] 


E 


0.00 


0.76 


1.49 


1.99 






A 


6.0 


0.25 


0.22 


1.86 




RPA [9, 10] 


A 


25.4 (22.8) 


<10 3 


0.04 (0.02) 


0.22 (0.74) 




Exp. [g, 26] 


E 


0.00 


0.95 


1.67 


2.62 






A 


6.00 ±0.40 


0.21 ±0.10 


0.18±0.10 


0.62 ±0.20 


38 ±1 



4. Final remarks 

We have shown that, when the capture of muons is evaluated in the context of the QRPA, the 
conservation of the number of particles is very important not only for carbon but in all light nuclei 
with A < 30. The consequence of this is the superiority of the PQRPA on the QRPA in this nuclear 
mass region, as can be seen from Fig. [lj 

The violation of the CVC by the Coulomb field in this mass region is of minor importance, 
since in (2.4) is AEVjoui ±-Eg is = 1 1.7 MeV, which is small in comparison with m^. However, this 
effect could be quite relevant for medium and heavy nuclei studied in Refs. [Q, O]. For instance, 
for 208 Pb is A£coui +Eg = 39.0 MeV, which implies a reduction of the operator 0@j — Oo,j for 
natural parity states by a factor 0.37, or equivalently that its contribution is only ~ 13% of that 
when the Coulomb field is not considered. 



We agree with the finding of Kortelainen and Suhonen Q22| ] on the extreme sensitivity of the 
muon capture rates on the pp coupling strength when described within the QRPA, as well as on 
a possible collapse of this approximation for the J^ = if state. Yet, in our opinion the QRPA 
behaves in this way dominantly in nuclei with a large neutron excess such as those analyzed in 



Refs. [ lj| , |13| ]. It 1S clear that the RQRPA calculation []T3p is sensitive to the pp coupling, while the 
RPA+BCS model [|l|] is not since it totally ignores the pp interaction. 

Finally, we conclude that the comparison between theory and data for the inclusive muon 
capture is not a fully satisfactory test on a nuclear model. The exclusive muon transitions are much 
more robust with respect to such a comparison. 
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